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Interplay between the dependence of symmetry energy on density and the variation of 
nucleonic densities across nuclear surface is discussed. That interplay gives rise to the 
mass dependence of the symmetry coefficient in an energy formula. Charge symmetry 
of the nuclear interactions allows to introduce isoscalar and isovector densities that are 
approximately independent of the magnitude of neutron-proton asymmetry. 

1. Introduction 

The utility of nuclear symmetry energy in assessing nuclear properties is a con- 
sequence of the charge symmetry of nuclear interactions, an invariance of those 
interactions under neutron-proton interchange. A broader symmetry is the charge 
invariance, an invariance under rotations in neutron-proton space. In the following, 
we shall explore different consequences of both types of invariance on nuclear char- 
acteristics, as well as interdependencies between those consequences. Mu ch o f the 
discussion is self-contained; additional technicalities can be found in RefsP^. 

One consequence of the charge symmetry is that isoscalar quantities can be 
identified within a nuclear system, that do not change under the proton-neutron 
interchange. These include the nuclear part of net energy and the net nucleon 
density, p(r) = p n (r)+p p (r). When expanded in relative asymmetry, rj — (N—Z)/A, 
an isoscalar quantity, F(rj), contains even powers of rj only, F(r]) = Fq + F2 rj 2 + 
F4T] 4 + . . . Due to the lack of a linear term, an isoscalar quantity depends weakly 
on asymmetry. Notably, for the expansion to apply, a nuclear quantity such as 
energy must be smoothed out to suppress shell and pairing effects. Apart from 
the isoscalar, isovector quantities may be identified, that change sign under the 
n <-> p interchange. An example is the neutron-proton density difference, p np (r) = 
Pn{ r ) —pp(r). When expanded in 77, an isovector quantity, G(rj), contains odd powers 
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of 77 only, G(j]) = G\ rj+G^ ?/ 3 +. . . Dividing an isovector quantity by 77, or by another 
isovector quantity, yields an isoscalar quantity, weakly dependent on 77: G/rj = 
Gi + G^ff + ■ ■ ■■ When considering the charge invariance, isovector quantities need 
to transform in a covariant fashion under a broader class of n-p transformations. 

Within the realm of consequences of charge-symmetry, we will be interested in 
the interplay of the density-dependence of symmetry energy in uniform matter and 
the characteristics of nucleonic densities across nuclear surface. We shall test our 
qualitative assertions within the Skyrme-Hartree-Fock (SHF) calculations of half- 
infinite nuclear matter that is the simplest nuclear system with a surface. In this 
context, we shall further consider mass-dependence of the symmetry coefficient in 
an energy formula. 



2. Symmetry Energy 

On account of the charge symmetry, the expansion of energy per nucleon of uniform 
matter, in asymmetry, takes the form 

j(Pn,P P ) = f (p) + S{p) ( B ^) 2 + ■ ■ (1) 

The leading term in this expansion is the energy per nucleon of symmetric matter. 
The remainder, the leading term of the remainder, and the coefficient S in the 
remainder are all called the symmetry energy in literature. Under most circum- 
stances, higher-order terms, beyond the first two on the r.h.s., play little role in 
the net energy. Correspondingly, the knowledge of two functions of density, ^ (p) 
and S(fl), allows, in practice, to describe the nuclear energy in a wide density and 
asymmetry range. E.g. in pure neutron matter, the nuclear energy is ^ ~ + S. 
The symmetry energy is customarily expanded around the normal density po, in a 
manner similar to the energy of symmetric matter, S(p) = a% + p ~^° + . . . Be- 
cause the energy of symmetric matter minimizes at po, the constant L determines 
the nuclear contribution to net pressure in the vicinity of p$, that is important for 
the structure of neutron stariM P ~ p 2 4^ ~ L . 

' 1 dp 3po 

In the Bethe-Weizsacker formula, expressing the energy of a nucleus in terms 
of nucleon numbers, the charge symmetry constraints the symmetry term to 
a quadratic form in asymmetry: 

Z 2 (N — Z) 2 
E = -a v A + a s A 2 / 3 + a c —+a a (A) i —L+E mic . (2) 

The symmetry coefficient a a can principally depend on mass number A, rather 
than having a fixed value {a a = a^) that is usually assumed. We shall explore 
a possibility of the mass dependence below. 

Due to its quadratic dependence on asymmetry, the nuclear contribution to the 
macroscopic energy of a nucleus, 

E = -a v A + a s A 2 ' 3 + ^ (N - Z) 2 = E (A) + ^ (N - Z) 2 , (3) 
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is analogous^ to the energy of a capacitor, E = Eq + Q 2 /2C, with the asymmetry 
N — Z being analogous to the charge Q and with the ratio A/2a a being analogous 
to the capacitance C . The asymmetry chemical potential, 

dE ^ 2a a (A) 

is, in this context, analogous to the electric potential V across capacitor terminals, 
V = Q/C. Turning back to uniform matter, we have 

_ dE _ 25(p) 
Ma ~ d(N -Z)~ Pnp p ' W 

and the ratio p/2S can be identified as capacitance per unit volume, cf. (jTJ) . In con- 
nection to the electrostatic analogy, we may note that, for connected independent 
capacitors, the net charge partitions itself in proportion to individual capacitances 
and the capacitances, otherwise, add up. 



3. Nuclear Densities 

We have mentioned that the net nucleon density p(r) = p n (r) + p p (r) is isoscalar 
and, as such, should weakly depend on (N — Z) for a given A. We will be ignoring 
here the Coulomb effects that may be accounted for in terms of corrections. The 
density difference, p np (r) = p n (f) — P P (r), is isovector but the renormalized density 
A p np (r) / (N — Z) is isoscalar in nature. Deficiency of the net relative asymmetry, 
as a renormalizing factor, is that it refers to global system characteristics. Quanti- 
tatively similar factor, that pertains to local properties only, is p, a /2a£ . With this, 
we introduce asymmetric density, a counterpart to the net density, as 

Pa(r) = ^[p n (r)-p P (r)] . (6) 

The asymmetric density serves as an isoscalar formfactor for the isovector density. 
With the specific normalization, the asymmetric density approaches normal density 
Po when the net density itself approaches po, in a weakly nonuniform system. 

As a consequence of charge symmetry, the nucleonic densities may be repre- 
sented in terms of two densities, p and p a , that each depends weakly on asymmetry^, 

Pn,p( r ) = \ [P( r ) ± ^?Pa( r )] ■ ( 7 ) 
L za a 

It is common to approximate the net nuclear density in terms of a Fermi shape 
p(r) = po/[l + cxp( I -j^)], where R = tqA 1 / 3 . Obviously, it is interesting to ask 
what the corresponding features of p a are. It turns out that those features are 
simultaneously related to a a (A) and S(p). 

Thus, the net capacitance for asymmetry may be written as 

A N-Z f Pnp 2 f 

= I dr — = — drp a (r). (8) 



2a a {A) p a J Pa 
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It follows then that the generalized asymmetry coefficient may be obtained by 
integrating the asymmetry density over volume. Upon subtracting and adding the 
isoscalar from the isovector density, we further find 

j£j - £ / dV(r) + ± J fir („. - ,)(,) K £ + ^ , (9) 

In the last step, we have noted that the two densities are close to po within the 
nuclear interior; thus, the last integration is concentrated to the surface region 
for a sufficiently large system, with the result being then proportional to A 2 / 3 . 
The surface modifies the capacity of a system for asymmetry, as compared to that 
expected for a uniform system at normal density. For large systems, the modification 
may be described in terms of the surface symmetry coefficient . 

Following ([5]), the neutron- proton density difference in a uniform system is 

Pnv = ^^ y (10) 

and the isovector density is then 

i v v 
_ 2a a p np _ a a p 

Pa — — a , n • I 11 ; 

Pa S{p) 

Given the short-range nature of nuclear interactions, these results should hold for 
weak nonuniformities, with the densities pertaining then to a specific location. In 
the context of (jTTJ) , we can rewrite the result ([8]) for the system capacitance as 

A ~ /" d 3 r J/W +.... (12) 



2a a {A) J 2S(p(r)) 
Here, the first r.h.s. term represents the regions where the local approximation (| 1 1 1) 
holds, and the dots represent the remainder. In (|12p . the net capacitance of the 
system emerges as an integral over independent capacitances at different localities; 
in pop the asymmetry is seen to partition itself in proportion to the local capaci- 
tance. In the context of (jlip and (d2), one can observe that dropping of symmetry 
energy in the surface region should act to enhance asymmetric density there and 
to increase the net system capacitance for asymmetry. The relative increase in net 
capacitance, cf. © , should be more pronounced in light systems where the surface 
plays a larger role than in heavy systems. 

4. Half-Infinite Matter in Skyrme-Hartree-Fock Calculations 

We next verify the validity of our qualitative claims by examining the features of 
half-infinite matter in the Skyrme-Hartree-Fock (SHF) calculations^^!. The direc- 
tion of uniformity is taken as z, with the vacuum and matter positioned at the 
positive and negative ends of the z-axis, respectively. The wavefunctions are taken 
in the form <f>(r) = 0(z)e 4ki ' ri , the wavevector space is discretized and we solve 
the equations for individual wavefunctions, 

- i £*u + (M) + u(z) ) m = e(k) Hz) ' (13) 
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where to* and U are obtained through self-consistency. 

We have claimed that the isoscalar and isovector densities should change little 
with asymmetry. The density profiles for the nuclear matter at different asymme- 
tries, for different interactions in SHF calculations, are shown in Figure [1] It is 
apparent that indeed either density changes little with r\. On the other hand, the 
isovector density changes significantly with interaction. 
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Fig. 1. Isocalar (left) and isovector (right) density profiles at different asymmetries in half-infinite 
nuclear matter from SHF calculations, for sample interactions considered in Refffl 



Figure [5] shows next the comparison of isoscalar and isovector densities for the 
different interactions. The larger the value of slope parameter L for the symmetry 
energy, the farther is the isovector density pushed out relative to the isoscalar 
density. This is consistent with the expectation based on the local approximation 
of Eq. The local approximation is also tested in Fig. [2] Up to the Friedel 

oscillations^, with a wavelength of A — ir/kp — 2.36 fm, the isovector density 
follows the local approximation down to the net density equal to about a quarter 
of normal density. An analysis^ following the WKB approximation shows that the 
local approximation may hold in the classically allowed region; in that respect, the 
quarter of normal density appears then to represent a typical classical return point. 

The farther the isovector density extends outside of the isoscalar density, with 
the larger L, the more significant is the modification of system capacitance by the 
surface, and the smaller is the surface symmetry coefficient af, cf. ([5]) and ©. 
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Fig. 2. Left panel shows the comparison of isoscalar and isovector densities, and of local approxima- 
tion to the isovector density, in symmetric half-infinite nuclear matter from the SHF calculations, 
for sample interactions. The longer- and shorter-dashed vertical lines, for a specific interaction, 
indicate, respectively, the location where the net density is equal to the quarter of normal density 
and the location of a classical return point for the Fermi wavevector directed along the z-axis. Right 
panel shows the ratio of symmetry energy coefficients a Y / , for different Skyrme interactions^, 
plotted vs the slope parameter L of the symmetry energy, scaled with the value of symmetry 
energy a%. 



The correlation between the ratios a^/af and L/a% is shown, for the Skyrme 
interactions, in the left panel of Fig. [2] Typical value of the volume coefficient for 
the interactions is aY ~ 30 MeV. 



5. Constraints from Isobaric Analog States 

We now turn to the constraints on symmetry energy stemming from an applica- 
tion of the energy formula (J2]). Unfortunately, in practice, competition between 
different physics terms within an energy formula makes it difficult to learn about 
the mass dependence of the symmetry coefficient, by fitting directly the formula 
to the ground-state nuclear energies^. However, the unwanted competition may be 
practically eliminated by exploiting the charge invariancc of nuclear interactions 
and generalizing the macroscopic energy formula to the lowest states of a given net 
isospin within a nucleus^- 1 . The generalization amounts to the replacement of the 
symmetry term in Eq. (|2|): 

a a (A) {N ~ Zf = 4«a(A) S 4a a (A) . {U) 
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In the ground state of a nucleus, the isospin takes on the lowest possible value T = 
\T Z \ — \N — Z\/2. This replacement absorbs the so-called Wigner term from E m i c . 
Following the formula generalization, it becomes possible to deduce the symmetry 
coefficient within a single nucleus, by using excitation energies to the isobaric analog 
states (IAS) representing the ground states of neighboring nuclei, with 

A(T(T + 1)) 

AE = 4a a (A)^^j '-± + AE mlc . (15) 

We use the tabulated energies of isobaric analog states^ and microscopic correc- 
tions to energies by Koura et aiP^to deduce the symmetry coefficients for individual 
nuclear masses. The mentioned corrections include deformation effects. Deduced co- 
efficient values range from a a ~ 9MeV, for light A < 10 nuclei, to a a ~ 22.5 MeV, 
for A > 200. Figure [3] shows inverse coefficient values plotted against inverse cube 
root of mass number. For A > 20, the value systematic is approximately linear, as 
expected from Eq. ((SJ). A fit with the r.h.s. of Eq. ([9]) produces coefficient values of 

= 32.9 MeV and of = 11.3 MeV. Similar coefficient values are obtained when 
trying to describe the a a (A)-data in terms of the Thomas- Fermi (TF) theorjI^ISl 
represented in the figure. The TF results suggest, nonetheless, that the effects of 
curvature of nuclear surface may, at some level, affect the symmetry energy even 
in the hea viest encountered nuclei. It should be mentioned that, within the TF 
theorjEEl, the local result II lip for the isovector density is exact. 




Fig. 3. Inverse of the symmetry coefficient, on the left scale, as a function of the inverse cube root 
of mass number, on the bottom scale. The filled circles represent the coefficients extracted from 
IAS excitation energies, with the microscopic corrections applied. The line shows a fit to the IAS 
results for A > 20. The squares represent the coefficients from a Thomas-Fermi model that best 
describes the data. 



8 Pawel Danielewicz and Jenny Lee 



4 — 



linear fit: 



3 



Thomas— Fermi 




03 



> <a 



2 — 



1 — 



Weizsacker 
mass fit 







O i 



20 



25 



30 



35 



40 



al (MeV) 



Fig. 4. Symmetry-energy parameter values in the plane of ajf /a„ and a£ , from comparing different 
theoretical expectations to data. 

The parameters arrived at in different ways are further summarized in Fig. [4] 
The favored values of symmetry energy at po are — (30 — 33) MeV. From Figs. |4] 
and [21 we read off the favored value of L ~ 80 MeV. 

Subsequent efforts in this investigation will be aimed at an extraction of the 
mass-dependent symmetry coefficients from spherical SHF calculations, at a model- 
independent understanding of the asymmetry skins, as well as at an understanding 
of the shell, Coulomb, deformation and curvature effects in the context of symmetry 
energy. 
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